Excitation and damping of collective modes of a Bose-Einstein condensate in a 

one-dimensional lattice 
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The mode structure of a Bose-Einstein condensate non- 
adiabatically loaded into a one-dimensional optical lattice is 
studied by analyzing the visibility of the interference pat- 
tern as well as the radial profile of the condensate after a 
time-of-flight. A simple model is proposed that predicts the 
short-time decrease of the visibility as a function of the con- 
densate parameters. In the radial direction, heavily damped 
oscillations are observed, as well as an increase in the con- 
densate temperature. These findings are interpreted as a re- 
thermalization due to dissipation of the initial condensate ex- 
citations into high-lying modes. 

PACS number(s): 03.75.Fi,32.80.Pj 



Studying collective modes in a Bose-Einstein conden- 
sate (BEC) is an efficient method for obtaining informa- 
tion about the dynamics of this quantum system 
So far, both experimental and theoretical results have 
been obtained for low-lying modes of a condensate in 
a magnetic trap, including breathing modes surface 
modes [Q| and the scissors mode ||^. Typically, in these 
experiments the collective modes were excited by a sud- 
den change in the trap frequency or geometry, and the 
frequency and damping rate of the subsequent oscilla- 
tions were measured either in situ or after a time-of-flight. 

In this paper we present experimental results and 
some preliminary theoretical considerations on the mode 
structure of a Bose-Einstein condensate inside a one- 
dimensional periodic potential. In experiments to date, 
BECs have been loaded into optical lattices mainly in 
the adiabatic regime in order to study, e.g., number 
squeezing and the Mott-insulator transition ||]. In 
this context, 'adiabatic' refers to the modes of the en- 
tire condensate rather than the single-well oscillation 
frequency. Therefore, the condition for adiabaticity is 
Tramp > whcrc Tramp IS the time over which 

the periodic potential is ramped up, and fi is the chem- 
ical potential of the condensate Q. By violating this 
condition (but still satisfying Tramp > '^T^/^^iat, with ujiat 
the single- well harmonic oscillator frequency), collective 
modes are excited in the condensate. 

Our exper imental setup is described in detail else- 
After creating BECs of A^o = 1-5(5) x 10* 
131 in a triaxial time-orbiting potential trap. 



trap an optical lattice along the vertical trap axis (for 
which the trap frequency viong ~ i^trap) created by two 
linearly polarized Gaussian laser beams intersecting at 
a half-angle 9 = 18 degrees and detuned by ~ 30 GHz 
above the rubidium resonance line. The periodic poten- 
tial V{z) = Vo sm^{nz/d) thus created has a lattice spac- 
ing d — 1.2 /im, and the depth Vq of the potential (mea- 
sured in lattice recoil energies Erec = fi^TT^ /2mcP) can 
be varied between and « 20 Erec by adjusting the laser 
intensity using acousto-optic modulators. The number of 
lattice sites occupied by the condensate lay between 10 
and 15, depending on the trap frequency. 




2 

P Prec 



where |10 
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we adiabatically lower the trap frequency Vtrap to the 
desired value and then superimpose onto the magnetic 
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FIG. 1. Evolution of the interference pattern (integrated 
perpendicular to the lattice direction) of a condensate re- 
leased from an optical lattice after non-adiabatic loading. 
The distinct two-peaked structure visible immediately after 
loading (a) is washed out within the first few milliseconds 
(b), after which the interference pattern takes on a com- 
plex structure ((c) and (d)). For long holding times, the ini- 
tial two-peaked structure reappears (e). In this experiment, 
T'trap = 26.7 Hz, Tramp = 5 ms, and Vo « 15Erec- From (a) 
to (e), thoid = 1,22,50, 100 and 300ms, respectively. 
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In a typical experiment, the optical lattice was ramped 
1 — 5 ms, after which the potential was 



up m Tj-djjip 

kept at its maximum value Vq for a holding time thoid- 
At the end of the holding time, the lattice was acceler- 
ated in 1 ms to one (lattice) recoil velocity by chirping the 
frequency difference between the lattice beams. Immedi- 
ately after that, both the lattice and the magnetic trap 
were switched off. After a time-of-fiight of 20 — 22 ms, 
the expanded condensate was imaged using a resonant 
probe flash. Figure ^ shows typical integrated absorp- 
tion images obtained in this manner for different holding 
times. For short times, a clean double-peak structure 
is visible, as expected from the interference between the 
condensates expanding from the individual lattice wells 
(with a TT phase difference between them due to the fi- 
nal acceleration). During the first few milliseconds, this 
pattern evolves into a more complicated structure fea- 
turing several additional peaks, and finally washes out 
completely, resulting in a single Gaussian-shaped lump. 
In an intermediate regime {thoid ~ 20 — 100 ms), the in- 
terference pattern again becomes complex, with multiple 
peaks. For longer waiting times, the two-peaked struc- 
ture reappears. 

We analyzed our experimental data in two different 
ways. In the lattice direction, we characterized the inter- 
ference pattern (integrated perpendicular to the lattice 
direction) through a visibility <^ calculated as 
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hpeak H~ ^middle 



(1) 



where hp^ak is the mean value of the absorption image 
at the position of the two peaks, and hmiddie is the value 
of the absorption image mid-way between the two peaks 
(averaged over a range of 1/5 of the peak separation). In 
the direction perpendicular to the optical lattice, we fit- 
ted a bimodal Gaussian function to the (longitudinally) 
integrated absorption profile and extracted from this fit 
the condensate fraction and the width pperp of the con- 
densate part. 
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FIG. 2. Examples of the short-term (a) and long-term be- 
haviour (b) of the visibility for a condensate in an optical 
lattice after non-adiabatic loading. In (a), the initial loss 
of visibility is evident (solid line: exponential fit with time 
constant 13ms). The trap frequency utrap = 19.8 Hz, and 
Vo/Ercc ~ 15 with Tramp = 1 ms. In (b), after strong ini- 
tial fluctuations the visibility approaches a stable value close 



to the initial visibility. In this experiment, l^tra 



26.7 Hz, 



Vo/Erec ~ 15 and Tramp = 5 ms. The dashed lines are an 
envelope to the data points with a constant upper part and 
an exponential lower part with a time constant of « 80 ms. 



Figure ^ shows examples of the short-term and long- 
term behaviour of ^. Initially, ^ rapidly decreases from 
a value of w 0.6 — 0.9 to roughly within « 5 — 20 ms, 
depending on the trap frequency vtrap- Subsequently, ^ 
typically rises again and begins to fluctuate in an ap- 
parently random manner for a few tens of milliseconds. 
Finally, these fluctuations die out and ^ stabilizes at a 
value close to the initial visibility. 

The corresponding behaviour of the condensate width 
Pperp and temperature (calculated from the condensate 
fraction) is shown in Fig. ||. One clearly sees radial os- 
cillations of Pperp at a frequency Vosc = '2i'perp, where 
V^i^trap is the trap frequency perpendicular 



to the lattice direction. These oscillations are heav- 
ily damped with a quality factor Q — ^irvosc/ldamp ~ 
10 — 30. Simultaneously, the condensate temperature ex- 
pressed as T/Tc (where Tj. is the critical temperature 
for the BEG transition) increases and approaches a new 
steady-state value on a time scale comparable to 1/ ^damp- 
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FIG. 3. Time evolution of the width (a) and the temper- 
ature (b) of the condensate in the optical lattice (same trap 
and lattice parameters as in Fig. ^ (b)). The solid lines are 
are fits to the data using an exponentially damped sinusoidal 
oscillation in (a) and a simple exponential function in (b). 
In both cases, the time constants of the exponential part are 

70 ms. Note the different time scales in the two graphs. 



The short-term decrease in ^ can be explained using 
a simple model. We take the Bose-Einstein condensate 
wavefunction to evolve according the Gross-Pitaevskii 
equation 
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Ui)sm2(^)+C/|^| 



(2) 



where Xj = {x, y, z} are the spatial coordinates along 
which the respective trap frequencies are Uj (= STTi^j), 
U = Anah"^ /m is the collisional interaction strength with 
a ~ 5.4 nm the s-wave scattering length. 

The rate at which the lattice is raised is slow enough 
that band excitations can be ignored and the condensate 
density will be reshaped to lie at the potential minima of 
each lattice site. Because we are interested in the phase 
properties of the entire condensate, it is convenient to 
consider a continuous Wannier or envelope represen- 
tation [T^ ] of the wavefunction, where we define the en- 
velope /(x, t) to represent the slowly varying amplitude 
of the Gross-Pitaevskii wavefunction V'(x, t), in which 
the rapid density variation along the lattice is smoothed 
out (also see M]). The evolution equation for / is of 
a similar form to Eq. ^ except that the lattice poten- 
tial no longer appears explicitly; the diffusion along z is 



generated by the Bloch dispersion relation reflecting the 
modified tunneling properties in this direction; and the 
mean- field term is renormalized as C/|'0P — * (l + c)C/|/p 
resulting from the compressed condensate density in the 
lattice. In the tight binding limit the Wannier states 
are localized and can be approximated by the harmonic 
oscillator orbital w{z) — exp(— z^/2a^g) / •(/ Tra^^, where 
aho = \/lT'lmujiat is the oscillator length and ujiat = 
•n^/^VQlmd? is the oscillation frequency about the lattice 
minima. In this limit the renormalized mean-field term 
is given by c= {d \\w{z)\'^ dz - X) = {{j mVod^ /2h^ - 1) 
Jirt , with the validity condition a <^ a^o ^ d. 

To develop an approximate expression for the short 
time dephasing behavior of the condensate we assume 
that during the lattice loading and subsequent time over 
which the dephasing occurs density transport in the con- 
densate is negligible, i.e. the evolution occurs in the 
phase of the wavefunction. In the lattice direction the 
tunneling frequency determines the time scale over which 
this assumption will be valid, which is typically of order 
100 ms in experiments where dephasing is observed. We 
also find that for a given trap frequency, there always ex- 
ists a minimum lattice depth below which no dephasing 
occurs. 

Expressing the envelope function in terms of amplitude 
and phase as / = |/| exp(iiS'), and neglecting the effects of 
spatial diffusion, the envelope equation of motion reduces 
to 



dt ^ ' 2 
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(7(1 + 5)1/1 = 



(3) 



Ignoring density transport |/p can be approximated 
by the Thomas-Fermi density profile for the conden- 
sate in the initial harmonic trap, which is of the form 
|V'tf(x)P = no[l — where no is the peak 

density and Rj = {R^, Ry, Rz} are the Thomas- Fermi 
radii With this substitution the phase evolves as 

where So(t) is a spatially constant phase. This result 
shows that the condensate develops a quadratic phase 
profile along the lattice arising from the imbalance of 
meanfield and harmonic trap potential energy in the lat- 
tice, and is equivalent to a spread in the quasi momentum 
distribution of the system. 
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FIG. 4. Momentum width of a condensate with a 
Thomas-Fermi density distribution and a quadratic phase 
profile. Numerical calculation (solid), linear fit (dotted). 

For a large lattice the quasi momentum distribution 
of the condensate and the momentum distribution of the 
envelope function are identical if the momentum distri- 
bution lies entirely within the first Brillouin zone. We 
have calculated the rms- width of the z-component of mo- 
mentum cTpz for a wavcfunction of the form ip{:x., (3) — 

■!/'tf(x) exp ^— i/3 (z/Rz)^^ , where /? is the total phase 

difference between the center and the outside of edge of 
the condensate. The results (shown in Fig. ^) indicate 
that for [3 2, the momentum width linearly increases 
with fi and is well approximated by apz{P) « 0.73 Ti(j/Rz. 
Identifying (3 with the coefficient of the quadratic spatial 
phase term in Eq. (^, we see that the width of the mo- 
mentum distribution increases linearly with time. 

The condensate will initially appear dephased when 
the quasi momentum width significantly fills the first 
Brillouin zone, which has a half width of h/2d. The ex- 
act portion of the Brillouin zone which must be filled de- 
pends on the observable used to determine the dephasing 
and may be taken as a fitting parameter. Here we take 
Cpz = h/Ad as the requirement for dephasing, which can 
be inverted using expression (Q) to yield the dephasing 
time 
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Figure ^ shows Tdeph as a function of viong (which deter- 
mines Rz and no) together with experimentally measured 
dephasing times. Equally good agreement with experi- 
ment is found when calculating Tdeph for the parameters 
of §. 
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FIG. 5. Experimental dephasing times as a function of the 
trap frequency viong in the lattice direction. The solid line is 
the theoretical prediction for Tdeph (see text). 

For the intermediate and long-term behaviour of the 
visibility, the phase-winding model predicts partial and 
complete revivals of the relative phases, leading to a com- 
plicated time-evolution of ^. Experimentally, we see a 
rather erratic behaviour of ^ for intermediate times, with 
the visibility fluctuation between — 0.2 and 0.5 — 0.6, 
followed by a stabilization at a high value around 0.6. 
The fact that this stabilization takes place on the same 
time-scale as the damping of the radial oscillations and 
the increase in temperature of the condensate leads us to 
speculate that all these phenomena are related through a 
dissipation mechanism whereby the the initial excitations 
in the longitudinal direction {i.e. reflected in the phase 
differences between adjacent lattice sites) and in the ra- 
dial directions lead to a re-thermalization of the system 
at a higher temperature and hence lower condensate frac- 
tion, with the remaining condensate now in the ground 
state of the combined harmonic and periodic trapping po- 
tentials. More detailed theoretical investigations in this 
direction are planned for the future. 

We also note here that the tunneling time between ad- 
jacent wells is also comparable to the damping time of 
the visibility fluctuations for the parameters of our ex- 
periment. In order to understand better the importance 
of the various mechanisms, i.e. damping with a possi- 
ble coupling between the different modes, and re-phasing 
due to tunneling, it will be important in future to conduct 
experiments for a range of different parameter combina- 
tions for the trap and the lattice. Also, repeating our 
experiment in a conflguration in which more lattice sites 
are occupied by the condensate will most likely eliminate 
the intermediate revivals (which require a large fraction 
of the sites to be in phase, which is less likely for a large 
number of sites) and leave only the lower envelope (see 
Fig. H) of the visibility evolution. Finally, the role of the 
flnite temperature (T « 0.7 Tc) at which we start our 
experiment will have to be investigated more closely. 

In summary, we have studied the behaviour of a Bose- 
Einstein condensate non-adiabatically loaded into a one- 
dimensional optical lattice. The time-evolution of the 
collective modes thus excited has been characterized 
through the visibility of the interference pattern as well 
as its radial width. Both quantities exhibit large ini- 
tial variations that are strongly damped on a time-scale 
comparable to the relaxation time of the condensate tem- 
perature. 
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